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ABSTRACT 

Anomalous Narrow-Line Quasars (ANLs) are a population of quasars with narrow 
H/3, and sometimes [Olll] broader than ~ 1000 km/s, in total comprising ~ 10 — 30% 
(most likely ~ 25%) of Type I quasars at 0.2 < z < 0.8. We hnd that virial masses 
using the H/3 and Mgll lines systematically differ for ANLs by an average of as much 
as 0.5 dex. Because the broad H/3 component width increases in ANLs but Mgll does 
not, we might suspect H/3-based virial masses for ANLs are wrong but Mgll masses 
are correct. If this is due to an outflow reaching the lower-ionization potential H/3 line, 
ClV masses will be similarly flawed. However, we cannot be certain of this explanation 
without followup work, and may be unable to identify which quasars are ANLs at 
z > 0.8. Therefore, it is essential that ANLs be well-understood and well-modeled in 
order to allow the use of virial mass estimators on large optical spectroscopic catalogs, 
particularly at z < 0.4 and z > 2.0 where only one broad line is available for use in 
mass estimation. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion disks 
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1 INTRODUCTION 

The quasar "broad-line region" (BLR) from which high- 
velocity gas produces correspondingly broad spectral emis- 
sion lines, with typical FWHM in the 2000-20000 km/s 
range. Prominent emission lines visible in optical spectra 
at z > 0.5 include, from ionization potential placing them 
nearest to the central black hole, CIV, a broad component 
of H/3, and Mgll. 

For some individual active galactic nuclei, rever bera- 
tion mapping (|Peterson et alj|2004l : iBentz et al.ll2009l ) has 
been able to confirm the locations of broad emission lines. 
In a time series of spectra for the same object, an increase 
in the continuum luminosity is followed, often hundreds of 
days later, by a similar flare in CIV and then H/3. Assum- 
ing the flare propagates outward at the speed of light, the 
delay can be used to infer a radius to BLR spectral lines. It 
is hoped that broad-line region velocities are predominantly 
virial. For reverberation-mapped quasars, this approxima- 
tion combined with Kepler's Laws then allows the use of 
CIV and H/3 to infer a black hole mass. 

An empirical relationship between the reverbera- 
tion radii a nd the continuum flux allows 'virial 
estimators' dMcLure fc Jarvk 
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iMcLure fc Dunlopl 
IWang et al.ri2009: 



Onken fc Kollmeierll2008l : iRisaliti et al.ll2009l : iRafiee fc Hani 
20101 ) that can be used on just one spectrum rather than 



a time series. Mass estimates have allowed a dramatic im- 
provement in our understanding of distant quasars in several 
ways, including the tight M — a correlation between black 
hole mass and stellar disp e rsion in the host gala ctic bulge 
l|Ferrarese fc Merrittl l200Ct iGebhardt etaH l200Ch , analysis 
of the quasar mass dist ribution as a function of redshift 
l|Vestergaard et all 20081), and analysis of the qua sar mass- 
luminosity plane (|Steinhardt fc Elvisll2010al . l201lh . The last 
two of these entirely depend upon virial mass estimates, 
while the others would require virial mass estimates if con- 
sidered at higher redshift. 

However, CIV in particular may also be substan- 
tially broadened by r adiation pressure and quasar outflows 
i|Marconi et al l |2009) . If so, this extra velocity would re- 
sult in an overestimate of the mass of the central black 
hole. It has recently been reported that ~ 25% of quasars 
at 0.2 < z < 0.8 have anomalous narrow lines, and that 
these ANL quasars are associated with broadened H/3 but 
not broadened Mgll. In this paper, we consider the implica- 
tions of this broadening on virial mass estimation. 

In § [2] we describe our line-fitting for the H/3/[OlIl] 
complex and identification of ANLs. We also examine the 
response of other emission lines to an increased narrow H/3 
width. The resulting response of virial mass estimates in 
ANLs is examined in § [3] Finally, the possible implications 
of these objects upon the validity of virial mass estimation 
are discussed in § [4] 
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Table 1. Properties of quasars in the SDSS DR7 catalog binned 
by narrow-line width 
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Figure 1. Comparison of the [OIII] FWHM and H/3 broad- 
component FWHM for quasars in the SDSS DR7 catalog. Lower 
narrow-line velocities are accompanied by a wide range of H/3 
widths, but anomalous narrow lines are well-correlated with broad 
H/3, increasing with increasing broad H/3 width. 



2 QUASARS WITH BROAD H/3 LINES 

W e investigate the propert i es of t he H/3 and [OIII] lines, as 
in ISteinhardt fc Silverman (l201lh . with a fitting prescrip- 
tion similar to IShen et al.l (|2008). First, the catalog red- 



shift is assumed and the continuum is fit at rest wave- 
lengths [4435, 4700] and [5100, 5535] A. The continuum is 
best-fit as a power law com bined with a Fe template from 
iBruhweiler fc Vernerl (^OOct ) , where the Fe template has two 
free parameters: (1) an amplitude and (2) convolution with 
a Gaussian of variable width. Once the continuum and iron 
lines are removed, H/3 line is fit with two Gaussians and 
the [OIII] doublet with a pair of Gaussians, such that the 
[OIII] doublet is constrained to have the 3:1 amplitude ratio 
physically required by the fine structure transitions involved 
and the two [OIII] widths are required to be identical to the 
width of the H/3 narrow component. 

A comparison between the H/3 broad-line FWHM and 
combined narrow H/3 and [OIII] FWHM shows two popula- 
tions (Fig. [T] a population with standard narrow lines un- 
corrected with the H/3 broad-line width and a population 
with broader "narrow" lines well correlated with the broad 
H/3 component. The key question for virial mass estimation, 
then, is whether this broader H/3 is due to an increased virial 
velocity or comes from non- virial broadening. We divide the 
SDSS DR7 sample by narrow- line width, as in Tabled We 
then co-add spectra within each bin to examine the depen- 
dence of Mgll emission on narrow H/3 width. Changes in 
other emission lines and the spectral continuu m of ANLs 
are discussed in ISteinhardt fc Silvermanl l |201lf ) but are be- 
yond the scope of this paper. 

Although the H/3 FWHM increases in ANLs, the Mgll 
line width (Fig. [2| does not increase. Since Mgll has an 
ionization potentia l placing it at a larger radius from the 
central black hole {Peterson 2008), one possible explanation 
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might be an outflow propagating only partway through the 
broad-line region. 

It is also important to determine whether other narrow 
lines might be good indicators of ANLs, since [OIII] is only 
visible at 0.2 < z < 0.8. For a the small fraction of our 
sample lying at z < 0.39, the [Sll] A6717, 6730 doublet is 
visible. In standard quasars, [Sll] is a pair of narrow lines, 
like [OIII] and the narrow component of H/3. Although very 
few objects are included at high [OIII] widths, it does appear 
that the [Sll] amplitude may be decreasing with increased 
narrow-line width, while the [Sll] width remains fixed (Fig. 
[2|. Still, [Sll] appears to be a poor indicator of whether a 
quasar is an ANL. 

However, [NeV] A3426 may be a better indicator, 
as it broadens with broad e r "na rrow" lines (Fig. [2j. As 
in ISteinhardt fc Silvermanl (|201ll ). approximately 11% of 
quasars are well-measured with broad [OIII] and a broad- 
ened H/3 narrow component, with an additional ~ 15% 
showing anomalous narrow lines where H/3 and [OIII] are 
poorly matched, most commonly with more broadening 
found in H/3 than in [OIII] . Thus, although [OIII] and [NeV] 
may turn out to be equally good indicators of ANLs, they 
are likely not good enough to identify the full set. It may, 
however, be that finding these two forbidden lines broadened 
together, while non-forbidden lines and one other forbidden 
line are not, is useful in developing a physical model for 
ANLs. 



3 ANLS AND VIRIAL MASS ESTIMATION 

One of the possible explanations for broader narrow H/3 
and [OIII], correlated with a broad H/3 component, is a 
strong outflow broadening all three lines. Such an outflow 
would result in substantial non-virial motion in the H/3 
broad component. Recent work using large catalogs to in- 
vestigate the properties of quasar ac cretion has been great!; 



enhanced by virial mass i n dicator s (McLurc & Jarvis 200 



iMcLure fc Dunlod 120041 : IVestergaard fc Peterson! 120061 1. 
predicated amongst other assumptions upon the idea that 
the gas producing H/3's broad component is virialized. 
This assumption has been questioned in more recent work 
|Wang et all 120091 : IRafiee fc Hallll2010l '), in which the virial 
velocity for broa d lines is fit a s v 2 = FWHM 7 rather than 
v 2 = FWHM 2 . IRafiee fc Hall (|2010l ) use various statisti- 
cal methods to find for Mgll a series of best-fit values 
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Figure 2. Comparison of co-added quasar spectra near the Mgll (top left), [SII] (top right), [OIII] (bottom left), and [NeV] (bottom 
right) emission lines at different narrow-line widths binned and colored as in Table[2] At each wavelength, only spectra at redshifts where 
those wavelengths can be measured were co-added. For [SII], there are too few objects at low z to produce a magenta curve. 



1.21 < 7 < 3.95, arguing that the best of these techniques 
is the one yielding 7 = 1.21 ± 0.40, within 2a of 7 = 2. 

A strong outflow would produce additional broaden- 
ing of H/3, and thus 7 < 2. Broader narrow lines are not, 
however, correlated with broadening in Mgll (Fig. [2]). Mgll 
and H/3 masses s how no offset, with a scatter of 0.22 dex 
ijShen et al.ll200Si ). Thus, if Mgll remains virial, a compar- 
ison between Mgll-based virial masses and H/3-based virial 
masses is a potential test for the nature of H/3 broaden- 
ing. Systematically larger H/3 masses in ANL systems would 
then indicate that H/3 is due to non-virial motion such as 
an outflow, while continued agreement would indicate that 
H/3 remains virial, perhaps broadened because ANL systems 
produce the proper ionization potential for H/3 at a different 
radius. 

The IShen et all (|201lh catalog pr oduces Mgll-based 
virial masses using the prescription of iMcLure fc Dunlopl 
l|2004h (along other estimators), which has been previously 
sho wn to agree with H /3 masses to within a scatter of ~ 0.22 
dex l|Shen et al.ll2008h . However, these catalogs restrict not 
just [OIII] but also the narrow component of H/3 to ^ 1200 
km/s in FWHM. ANLs include a broader H/3 narrow com- 



ponent, so restricting the FWHM to 1200 km/s will also 
underestimate the H/3 narrow component width. 

As a result, the line flux assumed by the fit to rep- 
resent the H/3 broad component would actually have been 
a combination of the broad component and the wings of 
the narrow component, resulting in a broad component fit 
that was artificially too narrow. An underestimate of the H/3 
FWHM will result in an underestimate in the virial mass 
Mbh oc FWHM 2 . 

Using the fitting routine described in § [2j we produce 
an alternative set of masses from the broad H/3 component, 
where the narrow component is no longer constrained to he 
at a FWHM less than 1200 km/s. For narrow [OIII] lines, 
there is good ag reement be t ween both sets of H/3 masses, but 
for ANLs, the lShen etail (|201ll ) masses are systematically 
lower (Fig. [3]). This i s likely due to the effect described above, 
where a IShen et al.l (|201lh fit a combination of the broad 
and narrow component as their broad component. Objects 
are only included in this sample if they have high enough 
signal-to-noise for our routine to produce a high-quality fit. 
Note that the high signal-to-noise sample contains a larger 
ANL fraction than the overall sample, since ANLs are on 
average more luminous than other quasars. 
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Figure 3. Comparison of H/3 masses between Shen et al. 2011, 
where the FWHM of [OIII] and narrow H/3 is constrained to be 
narrower than 1200 km/s, and masses produced in this paper, 
where these narrow lines are unconstrained. The sample consists 
of 10133 well-measured quasars at 0.4 < z < 0.8. Black points are 
objects with [OIII] narrower than 1200 km/s, while blue points 
correspond to objects with broader narrow components. 



Table 2. Mean virial mass comparison between this work and 
Shenll for well-measured quasars binned by [OIII] width, 0.5 < 
z < 0.6 
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For 4716 high signal-to-noise ANLs at 0.4 < z < 0.8 
(where both H/3 and Mgll masses are available), o ur H/3 
masses average 0.474 dex higher than IShen et al.l (|201ll ) 
using H/3, while for 5417 other quasars our masses aver- 
age 0.075 dex higher. Some of this discrepancy is due to 
changes in the Fell template used and other line-fitting de- 
tails, which are particula rly important for Mgll. Compared 
to the I Shen et al.l l|201ll ) Mgll masses, our masses average 
0.455 dex higher for ANLs and 0.170 dex higher for other 
quasars. Since the mass and luminosity d istributions change 
with redshift (|Steinhardt fc Elvisll2010al ). it is likely best to 
compare these mass estimates in a relatively narrow red- 
shift range, yet one wide enough to include a large enough 
sample to draw meaningful conclusions. At 0.5 < z < 0.6 
(Table [3j, our H/3 masses i|Vestergaard fc Peterson] |2006h 



average 0.047 dex higher than <|Shen et all 1201 ll ) for 1161 
"standard" quasars (FWHM < 1000 km/s), but 0.393 dex 
higher for 661 ANLs (> 1000 km/s). More generally, the 
discrepancy betw een our unconstrained H/3 masses and the 
IShen et all | 20Qg| ) Mgll masses, which are based upon a line 
that does not change in ANLs compared to other objects, 
grows with increasing [OIII] width (Table [3]). 

For comparis on, H/3 and Mgll virial masses in the 
Sh en et al.l (|2008l ) catalog show good agreement for all 
quasars. If the uncertainty in comparing H/3 and Mgll 
masses is uncorrelat ed between differe nt quasars with a scat- 
ter of ~ 0.22 dex ( Shen et al.l 12008? ) . the systematic error 
produced by restricting the narrow H/3 component to 1200 
km/s would be the dominant source of error in ANL masses. 

Alternative estimate s for the statistical uncertain ty in 
H/3 alone include 0.4 dex (IVestergaard fc Petersonll2006l ) and 
0.21 dex (|Steinhardt fc Elvis! l2010bl ). with the latter in ad- 
dition to 0.15 dex scatter in Mgll. Neither of these esti- 
mates investigated whether errors in H/3 and Mgll masses 
are correlated, as they will be for errors emanating from 
the empirical relationship between monochromatic contin- 
uum luminosity and inferred radius to the broad-line region. 
Even if the uncertainty were 0.4 dex and entirely statistical, 
the discrepancy between H/3 and Mgll mass estimates for 
ANLs would be at least comparable. Thus, this systematic 
disagreement cannot be neglected for individual objects, let 
alone in the analysis of large catalogs where the sample size 
allows a substantial reducti on in the statistic al error. 

We also note that the IShen et al.l l|201ll ) masses show 
much better agreement between H/3 Mgll for ANLs, despite 
requiring as part of the fitting routine that the narrow com- 
ponent of FL/3 along with the [OIII] doublet have FWHM 
^ 1200 km/s. As above, an artificially-low broad H/3 width 
will produce an artificially low virial mass M V - 1T oc FWHM 2 . 
Remarkably, this lower H/3 mass turns out to be empirically 
almost exactly the a mount needed to bring H/3 and Mgll 
back into agreement: IShen et al.l (|201ll ) H/3 masses average 
0.094 dex higher for 5417 well-measured "standard" quasars 
(FWHM < 1000 km/s), but 0.075 dex lower for 4716 well- 
measured ANLs ( > 1000 km/s). Similarly, at 0.5 < z < 0.6, 
IShen et al.l (|201ll ) H/3 masses average 0.111 dex higher for 
1161 "standard" quasars (FWHM < 1000 km/s), but 0.065 
dex lower for 661 ANLs (> 1000 km/s). It should be noted 
that although H/3 and Mgll masses appear close to agree- 
ment for both standard and ANL quasars, the typical offset 
has changed by ~ 0.17 dex between the two samples. 

Regardless, these are not based upon the full FWHM 
of the H/3 line. As demonstrated by iRafiee fc Halll l|2010l ). 
without a theoretical basis for choosing a mass estimator, 
different empirical calibration techniques will produce a wide 
variety of possible masses due to the very small number of 
reverberation mapping systems available. There is a theo- 
retical basis for using the entire line FWHM and assuming 
virial motion in order to determine the virial component of 
the broad-line region velocity as v v i r oc FWHM 2 . However, 
this relation now seems difficult to reconcile with the exis- 
tence of ANLs even for H/3, let alone CIV at a smaller radius 
from the central black hole. It may be possible to develop 
a theoretical basis for reducing the measured FWHM to a 
purely virial FWHM if the broadening can be well-modeled. 
However, we currently l ack that understa nding, and there 
is no basis for using the IShen et al.l (|201 lh prescription for 
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reducing the FWHM over any other that could be produced 
from calibrating H/3 masses against Mgll for ANLs. 

Virial mass estimation yields systematically different re- 
sults for ANLs when using H/3 and Mgll, even th ough the ef- 
fect w as previously masked by constraints in the lShen et al.l 
(2008) line fitting routine. At least one of the two must be 
wrong, either because of substantial non- virial motion in the 
broad line or because the empirical relationship between the 
radius to the broad-line region and continuum luminosity 
breaks down for ANLs. Because the Mgll line does not ap- 
pear to change with increasing [OIII] while the broad com- 
ponent of H/3 is strongly correlated, we might suspect that 
the H/3 masses are in error, perhaps due to a strong outflow. 
However, it is proper to be cautious regarding Mgll masses 
as well until there are better measurements, ideally includ- 
ing reverberation mapping, on ANLs and these calibrations 
can be examined more closely. 



4 DISCUSSION 

Virial mass estimates for anomalous narrow-line quasars 
(ANLs), a new population comprising perhaps one quarter 
of Type I quasars at 0.2 < z < 0.8, appear to be flawed. 
For ANLs, the assumption that the line FWHM is entirely 
due to the virial velocity of the gas producing H/3 and Mgll 
leads to a systematic disagreement in estimated mass, with 
the disagreement intensifying with broader [OIII]. We must 
conclude that at least one of these estimates is wrong. 

The prognosis for virial mass estimation is potentially 
dire. The best-case scenario is likely that H/3 is indeed broad- 
ened by an outflow, and we might still hope that Mgll, which 
appears unresponsive to that outflow, remains virial and 
can be used to produce a reliable mass. Note, however, that 
ANLs represent ~ 30% of quasars at 0.2 < z < 0.8. There- 
fore, at 0.2 < 2 < 0.4, where Mgll does not lie within the 
SDSS spectrograph, 30% of quasars will have an incorrect 
mass. At present, that mass will be uncorrectable, although 
a sufficiently precise model for the outflow may allow a fu- 
ture correction. 

It is at least fortunate that the presence of [OIII] (and 
the narrow H/3 component) around the H/3 line provides 
an indicator for many (but not all) of the 30% of objects 
that will have incorrect masses. However, much of the util- 
ity of virial masses has come from our ability to take a large, 
statistically-complete sample (such as provided by SDSS at 
i < 19.2) and look at the mass distribution, often compared 
against other quasar properties. ANLs are not randomly- 
selected quasars, but rather quasars with specific emission 
line properties, continuum properties, and even host galac- 
tic properties. Even removing all ANLs from the sample and 
only using the remaining 70% of quasars will likely result in 
incorrect conclusions. 

At high redshift , the problem bec omes even worse. Al- 
though the BOSS (jRoss et all boiOT l survey will provide 
a large, statistically-complete sample at high redshift, for 
z > 2.0, Mgll is no longer visible and CIV masses must 
be used. CIV has an ionization potential placing it even 
closer to the central black hole than H/3 and is more vulner- 
able to effects such a s radiation pressure and quasar wind 
ijMarconi et al J l2009h . and thus any outflow affecting H/3 
very likely will also affect CIV. So, we might expect ^ 30% 



of CIV masses to also be wrong due to outflows. Yet, for 
H/3, the presence of the [OIII] narrow line at least allows a 
detection of ANLs and confidence that non- ANLs will have 
better-measured masses. Although other narrow lines ex- 
ist in quasar spectra, [Sll] is not broadened in ANLs, even 
though [NeV] is. It is plausible that all narrow lines visible 
in the optical at z > 2.0 will remain narrow, with the strong 
semi-forbidden line CIII] A1909 perhaps our best hope for 
an ANL indicator. If so, not only will ^ 30% of CIV masses 
be overestimates, but we may also have no indication as 
to which masses are erroneous. Understanding ANLs well 
enough to model and outflows and correct for them would 
then become a requirement in order to make use of high- 
redshift quasar masses. 

The authors would like to thank Steve Balbus, Tim 
Brandt, Forrest Collman, Martin Elvis, Eilat Glikman, 
Jeremy Goodman, Julian Krolik, Greg Novak, Jerry Os- 
triker, David Spergel, Michael Strauss, and Todd Thompson 
for valuable comments. This work was supported by World 
Premier International Research Center Initiative (WPI Ini- 
tiative), MEXT, Japan. 



REFERENCES 

Bentz M. C, Peterson B. M., Netzer H., Pogge R. W., 

Vestergaard M., 2009, ApJ, 697, 160 
Bruhweiler F., Verner E., 2008, ApJ, 675, 83 
Ferrarese L., Merritt D, 2000, ApJ, 539, L9 
Gebhardt K., Kormendy J., Ho L. C, et al., 2000, ApJ, 

543, L5 

Marconi A., Axon D. J., Maiolino R., et al., 2009, ApJ, 
698, L103 

McLure R. J., Dunlop J. S., 2004, MNRAS, 352, 1390 
McLure R. J., Jarvis M. J., 2002, MNRAS, 337, 109 
Onken C. A., Kollmeier J. A., 2008, ApJ, 689, L13 
Peterson B., 2008, An Introduction to Active Galactic Nu- 
clei, Cambridge University Press: Cambridge 
Peterson B. M., Ferrarese L., Gilbert K. M., et al., 2004, 

ApJ, 613, 682 
Ranee A., Hall P. B., 2010, ArXiv:1011.1268 
Risaliti G., Young M., Elvis M., 2009, ApJ, 700, L6 
Ross N., Sheldon E. S., Myers A. D, et al., 2010, in Bul- 
letin of the American Astronomical Society, vol. 41 of 
Bulletin of the American Astronomical Society, 517 
Shen Y., Greene J. E., Strauss M. A., Richards G. T., 

Schneider D. P., 2008, ApJ, 680, 169 
Shen Y., Richards G. T., Strauss M. A., et al., 2011, ApJS, 
194, 45 

Steinhardt C. L., Elvis M., 2010a, MNRAS, 402, 2637 
Steinhardt C. L., Elvis M., 2010b, MNRAS, 406, LI 
Steinhardt C. L., Elvis M., 2011, MNRAS, 410, 201 
Steinhardt C. L., Silverman J. D, 2011, ArXiv:1109.0537 
Vestergaard M., Fan X., Tremonti C. A., Osmer P. S., 

Richards G. T., 2008, ApJ, 674, LI 
Vestergaard M., Peterson B. M., 2006, ApJ, 641, 689 
Wang J., Dong X., Wang T., et al., 2009, ApJ, 707, 1334 



